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Abstract 



We provide an up to date account of progress in understanding the origin 
of the observed baryon asymmetry of the universe. While our primary goal 
is to be current, we have attempted to give a pedagogical introduction to 
CN ' the primary areas of research in this field, giving a detailed description of 

■ the different scenarios. The very recent developments in GUT baryogenesis, 

■ leptogenesis, electroweak baryogenesis and the Affleck-Dine mechanism are 
presented. In particular, we focus on specific particle physics implementa- 



| tions, mostly in the context of supersymmetry, which lead to specific testable 

predictions. 

^ I. INTRODUCTION 

The symmetry between particles and antiparticles |],||]], firmly established in collider 
physics, naturally leads to the question of why the observed universe is composed almost 
entirely of matter with little or no primordial antimatter. 

Outside of particle accelerators, antimatter can be seen in cosmic rays in the form of a 
few antiprotons, present at a level of around 10~ 4 in comparison with the number of protons 
(for example see 0). However, this proportion is consistent with secondary antiproton 
production through accelerator-like processes, p + p — > 3p + p, as the cosmic rays stream 
towards us. Thus there is no evidence for primordial antimatter in our galaxy. Also, if matter 
and antimatter galaxies were to coexist in clusters of galaxies, then we would expect there 
to be a detectable background of 7-radiation from nucleon-antinucleon annihilations within 
the clusters. This background is not observed and so we conclude that there is negligible 
antimatter on the scale of clusters (For a review of the evidence for a baryon asymmetry see 
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More generally, if large domains of matter and antimatter exist, then annihilations would 
take place at the interfaces between them. If the typical size of such a domain was small 
enough, then the energy released by these annihilations would result in a diffuse 7-ray 
background and a distortion of the cosmic microwave radiation, neither of which is observed. 
A careful numerical analysis of this problem || demonstrates that the universe must consist 
entirely of either matter or antimatter on all scales up to the Hubble size. It therefore seems 
that the universe is fundamentally matter-antimatter asymmetric. 

While the above considerations put an experimental upper bound on the amount of an- 
timatter in the universe, strict quantitative estimates of the relative abundances of baryonic 
matter and antimatter may also be obtained from the standard cosmology. The baryon 
number density does not remain constant during the evolution of the universe, instead scal- 
ing like a -3 , where a is the cosmological scale factor ||. It is therefore convenient to define 
the baryon asymmetry of the universe in terms of the quantity 

V = — , (1-1) 

n 7 

where % = rif, — ni is the difference between the number of baryons and antibaryons per 
unit volume and n 7 = is the photon number density at temperature T. Primordial 

nucleosynthesis (for a review see 0) is one of the most powerful predictions of the standard 
cosmological model. The theory allows accurate predictions of the cosmological abundances 
of all the light elements, H, 3 He, 4 He, D, B and 7 Li, while requiring only the single input 
parameter 77 which has been constant since nucleosynthesis. The range of 77 consistent with 
the deuterium and 3 He primordial abundances is || 

4(3) x 1(T 10 < r] < 7(10) x 1(T 10 , (1.2) 

where the most conservative bounds are in parentheses. Alternatively we may write the 
range as 

0.015(0.011) < tt B h 2 < 0.026(0.038) , (1.3) 

where Qb is the proportion of the critical energy density in baryons, and 0.5 < h < 0.9 
parametrizes the present value of the Hubble parameter via h = H / (100 Km Mpc -1 sec -1 ). 

To see that the standard cosmological model cannot explain the observed value of 77, 
suppose that initially we start with rj = 0. We can compute the final number density of 
nucleons b that are left over after annihilations have frozen out. At temperatures T <> 1 
GeV the equilibrium abundance of nucleons and antinucleons is |J 

rth nr. /m„\ 3 / 2 m v 

-? 1.4 
n 7 n 7 \ J / 

When the universe cools off, the number of nucleons and antinucleons decreases as long 
as the annihilation rate r ann ~ ^(a^v) is larger than the expansion rate of the universe 
H ~ 1.66 gl^ 2 -^-. The thermally averaged annihilation cross section (o^f) is of the order of 
m 2 , so at T ~ 20 MeV, r ann ~ H, and annihilations freeze out, nucleons and antinucleons 
being so rare that they cannot annihilate any longer. Therefore, from ( |1.4| ) we obtain 
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~ ~ 10 , (1.5) 

which is much smaller than the value required by nucleosynthesis. 

In conclusion, in the standard cosmological model there is no explanation for the value 
of the ratio ( |1.2| ), if we start from rj = 0. An initial asymmetry may be imposed by hand as 
an initial condition, but this would violate any naturalness principle. Rather, the guiding 
principle behind modern cosmology is to attempt to explain the initial conditions required 
by the standard cosmology on the basis of quantum field theories of elementary particles in 
the early universe. In this context, the generation of the observed value of 7] is referred to 
as baryogenesis. 

The goal of this article is to describe recent progress in understanding scenarios for 
generating the baryon asymmetry of the universe (BAU), as quantified in equation fll.2|) , 
within the context of modern cosmology. Earlier reviews of this subject, with a different 
focus, can be found in Refs. f8HT3fl. We provide a brief account of the central ideas of 



baryogenesis, and then review the most recent views on the major implementations. We 
hope that this provides an up to date summary of the status of the field. 

In the next section we shall describe the Sakharov Criteria which must be satisfied by any 
particle physics theory through which a baryon asymmetry is produced. In section (|111| ) we 
describe how bayogenesis can proceed in certain Grand Unified theories and, in particular, 
we focus on how new ideas concerning preheating after inflation have modified this scenario. 
In section ( [TVD we review electroweak baryogenesis, concentrating on the implications of 
new results concerning the order of the electroweak phase transition (EWPT), the Higgs 
mass, and supersymmetric implementations. Section (0) contains an account of Affleck- 
Dine (AD) type baryogenesis scenarios. After describing the general ideas, we explain how 
this mechanism can find a natural home in the context of supersymmetric theories. We 
conclude by discussing the current status of baryogenesis and prospects for future progress. 

A note about our conventions. Throughout we use a metric with signature +2 and, 
unless explicitly stated otherwise, we employ units such that h — c — k — 1. 



II. THE SAKHAROV CRITERIA 



As pointed out by Sakharov [14| , a small baryon asymmetry rj may have been produced in 
the early universe if three necessary conditions are satisfied: i) baryon number (B) violation; 
ii) violation of C (charge conjugation symmetry) and CP (the composition of parity and 
C) and Hi) departure from thermal equilibrium. The first condition should be clear since, 
starting from a baryon symmetric universe with r] = 0, baryon number violation must take 
place in order to evolve into a universe in which r] does not vanish. The second Sakharov 
criterion is required because, if C and CP are exact symmetries, then one can prove that 
the total rate for any process which produces an excess of baryons is equal to the rate of 
the complementary process which produces an excess of antibaryons and so no net baryon 
number can be created. That is to say that the thermal average of the baryon number 
operator B, which is odd under both C and CP, is zero unless those discrete symmetries 
are violated. CP violation is present either if there are complex phases in the lagrangian 
which cannot be reabsorbed by field redifinitions (explicit breaking) or if some Higgs scalar 
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field acquires a VEV which is not real (spontaneous breaking). We will discuss this in detail 
shortly. 

Finally, to explain the third criterion, one can calculate the equilibrium average of B 

(B) T = Tr (e~ pH B) = Tr [{CPT)(CPT)- x e- pH B)\ 

= Tr (e-^iCPT^BiCPT)] = -Tr (e~ pH B) , (2.1) 

where we have used that the Hamiltonian H commutes with CPT. Thus (B)t = in 
equilibrium and there is no generation of net baryon number. 

Of the three Sakharov conditions, baryon number violation and C and CP violation may 
be investigated only within a given particle physics model, while the third condition - the 
departure from thermal equilibrium - may be discussed in a more general way, as we shall 
see. Let us discuss the Sakharov criteria in more detail. 



A. Baryon Number Violation 

1. B -violation in Grand Unified Theories 



Grand Unified Theories (GUTs) |15| describe the fundamental interactions by means 
of a unique gauge group G which contains the Standard Model (SM) gauge group 
SU(3)c ® SU(2)l <8> U(1)y- The fundamental idea of GUTs is that at energies higher than 
a certain energy threshold M GU t the group symmetry is G and that, at lower energies, the 
symmetry is broken down to the SM gauge symmetry, possibly through a chain of symmetry 
breakings. The main motivation for this scenario is that, at least in supersymmetric models, 
the (running) gauge couplings of the SM unify ||T6"|-|T8|1 at the scale Mgut — 2 x 10 16 GeV, 



hinting at the presence of a GUT involving a higher symmetry with a single gauge coupling. 

Baryon number violation seems very natural in GUTs. Indeed, a general property of 
these theories is that the same representation of G may contain both quarks and leptons, 
and therefore it is possible for scalar and gauge bosons to mediate gauge interactions among 
fermions having different baryon number. However, this alone is not sufficient to conclude 
that baryon number is automatically violated in GUTs, since in some circumstances it is 
possible to assign a baryonic charge to the gauge bosons in such a way that at each boson- 
fermion-fermion vertex the baryon number is conserved. In the particular case of the gauge 
group 377(5), it turns out that among all the scalar and gauge bosons which couple only to 
the fermions of the SM, five of them may give rise to interactions which violate the baryon 
number. The fermionic content of SU(5) is the same as that of the SM. Fermions are 
assigned to the reducible representation 5/ © 10/ as 5/ = {d c L ,£ L } and 10/ = {QL,u c L ,e c L }. 
In addition, there are 24 gauge bosons which belong to the adjoint representation 24y and 
which may couple to the fermions through the couplings 

(g/V2)24 v [(B,)t (E f ) + (10,)* 10/] . (2.2) 

Because of this structure, the gauge bosons XY with SM gauge numbers (3,2, —5/6) may 
undergo baryon number violating decays. While in the gauge sector of SU(5) the structure 
is uniquely determined by the gauge group, in the Higgs sector the results depend upon 
the choice of the representation. The Higgs fields which couple to the fermions may be in 
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the representation 5# or in the representations 10h, 15#, 45# and 50#. If we consider the 
minimal choice 5#, we obtain 

h n (10/f (10/) 5 H + h D (5/) T (10/) 5 H , (2.3) 

where Hjj,d are complex matrices in the flavor space. The representation 5# contains the 
Higgs doublet of the SM, (1,2,1/2) and the triplet H 3 = (3,1,-1/3) which is B-violating 
Note that each of these bosons have the same value of the charge combination B — L, 



which means that this symmetry may not be violated in any vertex of SU(5). 

However, if we extend the fermionic content of the theory beyond that of minimal SU(5), 
we allow the presence of more heavy bosons which may violate B and even B — L. In GUTs 
based on SO (10), for instance, there exists a fermion which is a singlet under the SM gauge 
group, carries lepton number L — — 1, and is identified with the right-handed antineutrino. 
However, this choice for the lepton number assignment leads to no new gauge boson which 
violates B — L. As we mentioned, the generation of a baryon asymmetry requires C violation. 
Since 50(10) is C-symmetric, the C symmetry must be broken before a baryon asymmetry 
may be created. If 50(10) breaks down to SU(2) L <g> SU(2) R <g> SU(4) via the nonzero 
expectation value of a Higgs field in the 54 H representation p0|-p2|, the C symmetry is not 
broken until U(1)b-l is broken at the scale Mb-l- In this case, the right-handed neutrino 
acquires a Majorana mass Mn = O(Mb-l) and its out-of-equilibrium decays may generate 
a nonvanishing B — L asymmetry [^3j . 



B. S-violation in the Electroweak theory. 

It is well-known that the most general Lagrangian invariant under the SM gauge group 
and containing only color singlet Higgs fields is automatically invariant under global abelian 
symmetries which may be identified with the baryonic and leptonic symmetries. These, 
therefore, are accidental symmetries and as a result it is not possible to violate B and 
L at tree-level or at any order of perturbation theory. Nevertheless, in many cases the 
perturbative expansion does not describe all the dynamics of the theory and, indeed, in 1976 
't Hooft |E3 realized that nonperturbative effects (instantons) may give rise to processes 



which violate the combination B + L, but not the orthogonal combination B — L. The 
probability of these processes occurring today is exponentially suppressed and probably 
irrelevant. However, in more extreme situations - like the primordial universe at very high 
temperatures |p5|- p^1 - baryon and lepton number violating processes may be fast enough 



to play a significant role in baryogenesis. Let us have a closer look. 



At the quantum level, the baryon and the lepton symmetries are anomalous |29| , |30 

f 

V32tt 2 "^" 32tt 



d,f B = d,f L = n, (j^W« u W^ - 4^ F ^ V ) - ( 2 - 4 ) 



where g and g' are the gauge couplings of SU(2) L and U(l)y, respectively, n/ is the number 
of families and = (l/2)e^ ua(3 W a(3 js the dual of the SU(2) L field strength tensor, with 
an analogous expression holding for F. To understand how the anomaly is closely related 
to the vacuum structure of the theory, we may compute the change in baryon number from 
time t = to some arbitrary final time t — tf. For transitions between vacua, the average 
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values of the field strengths are zero at the beginning and the end of the evolution. The 
change in baryon number may be written as 



AB = AN CS = n f [N cs (t f ) - N cs (0)] . (2.5) 
where the Cherns-Simon number is defined to be 

N CS (t) = 3^3 / d 3 xe^ k Tr [a^Aj, + hgA^A^ . (2.6) 

Although the Chern-Simons number is not gauge invariant, the change ANcs is- Thus, 
changes in Chern-Simons number result in changes in baryon number which are integral 
multiples of the number of families rif. Gauge transformations U(x) which connects two 
degenerate vacua of the gauge theory may change the Chern-Simons number by an integer 
n, the winding number. If the system is able to perform a transition from the vacuum Q^. to 
the closest one Q^t l \ the Chern-Simons number is changed by unity and AB = AL — rif. 
Each transition creates 9 left-handed quarks (3 color states for each generation) and 3 left- 
handed leptons (one per generation). However, adjacent vacua of the electroweak theory are 
separated by a ridge of configurations with energies larger than that of the vacuum. The 
lowest energy point on this ridge is a saddle point solution to the equations of motion with 
a single negative eigenvalue, and is referred to as the sphaleron p6|p7| . The probability of 
baryon number nonconserving processes at zero temperature has been computed by 't Hooft 
|2~3| and is highly suppressed by a factor exp(— 4n/a w ), where aw = g 2 /4ir. This factor 



may be interpreted as the probability of making a transition from one classical vacuum to 
the closest one by tunneling through an energy barrier of height ~ 10 TeV corresponding 
to the sphaleron. On the other hand, one might think that fast baryon number violating 
transitions may be obtained in physical situations which involve a large number of fields. 
Since the sphaleron may be produced by collective and coherent excitations containing ~ 
1/ctw quanta with wavelength of the order of 1/Myy, one expects that at temperatures T 3> 
Mw, these modes essentially obey statistical mechanics and the transition probability may 
be computed via classical considerations. The general framework for evaluating the thermal 
rate of anomalous processes in the electroweak theory was developed in |[28|| . Analogously 
to the case of zero temperature and since the transition which violates the baryon number 
is sustained by the sphaleron configuration, the thermal rate of baryon number violation in 
the broken phase is proportional to exp(— S3/T), where S3 is the three-diemensional action 
computed along the sphaleron configuration, 

S 3 = E sp (T) ee (M w (T)/a w )S , (2.7) 

with the dimensionless parameter £ lying in the range 3.1 < £ < 5.4 depending on the Higgs 
mass. The prefactor of the thermal rate was computed in |yj as 

r w (T)=^(^|) 3 Miexp(-^)) , (2.8) 

where \l is a dimensionless constant. Recent approaches to calculating the rate of baryon 
number violating events in the broken phase have been primarily numerical. The best 



calculation to date of the broken phase sphaleron rate is undoubtably that by Moore [32 
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This work yields a fully nonperturbative evaluation of the broken phase rate by using a 
combination of multicanonical weighting and real time techniques. 

Although the Boltzmann suppression in ( |2.8| ) appears large, it is to be expected that, 
when the electroweak symmetry becomes restored at a temperature of around 100 GeV, 
there will no longer be an exponential suppression factor. Although calculation of the 
baryon number violating rate in the high temperature unbroken phase is extremely difficult, 
a simple estimate is possible. The only important scale in the symmetric phase is the 
magnetic screening length given by £ = (awT)' 1 . Thus, on dimensional grounds, we expect 
the rate per unit volume of sphaleron events to be 

T sp (T) = n(a w T) 4 , (2.9) 

with k another dimensionless constant. The rate of sphaleron processes can be related to 



the diffusion constant for Chern-Simons number by a fluctuation-dissipation theorem |SB| 



(for a good description of this see [10]). In almost all numerical calculations of the sphaleron 
rate, this relationship is used and what is actually evaluated is the diffusion constant. The 
first attempts to numerically estimate k in this way @,[3lJ yielded k ~ 0.1 — 1, but the 



approach suffered from limited statistics and large volume systematic errors. Nevertheless, 
more recent numerical attempts p6|-[39|[ found approximately the same result. However, 



these approaches employ a poor definition of the Chern-Simons number which compromises 
their reliability. 



This simple scaling argument leading to (|2.9|) has been recently criticized in refs. [40,41 



where it has been argued that damping effects in the plasma suppress the rate by an extra 
power of aw to give Y sp ~ a w T 4 . Indeed, since the transition rate involves physics at 
soft energies g 2 T that are small compared to the typical hard energies ~ T of the thermal 
excitations in the plasma, the simplest way of analyzing the problem is to consider an 
effective theory for the soft modes, where the hard modes have been integrated out, and 
to keep the dominant contributions, the so-called hard thermal loops. It is the resulting 
typical frequency u c of a gauge field configuration immersed in the plasma and with spatial 
extent (g 2 T)~ l that determines the change of baryon number per unit time and unit volume. 
This frequency u c has been estimated to be ~ g A T when taking into account the damping 
effects of the hard modes f4"0| , |4"I |. The effective dynamics of soft nonabelian gauge fields 



at finite temperature has also been recently addressed in f42|j43|| , where it was found that 



Tsp ~ ce w T \n(l/a\y). Lattice simulations with hard-thermal loops included have been 
performed [Q and seem to indicate the T sp ~ 30a w T 4 , which is not far from ct^T 4 . 



1. Other ways of achieving B -violation. 



One of the features which distinguishes supersymmetric field theories from ordi- 

nary ones is the existence of "flat directions" in field space on which the scalar potential 
vanishes. At the level of renormalizable terms, such flat directions are generic. Supersymme- 
try breaking lifts these directions and sets the scale for their potential. ^From a cosmological 
perspective, these flat directions can have profound consequences. The parameters which 
describe the flat directions can be thought of as expectation values of massless chiral fields, 
the moduli. Many flat directions present in the minimal supersymmetric standard model 
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(MSSM), carry baryon or lepton number, since along them squarks and sleptons have non 
zero vacuum expectation values (VEVs). If baryon and lepton number are explicitly broken 
it is possible to excite a non-zero baryon or lepton number along such directions, as first 
suggested by Affleck and Dine |47j] in a more general context. We shall explore these ideas 
more thoroughly in section |V|. 

C. CP violation 

1. CP-Violation in Grand Unified Theories. 

CP violation in GUTs arises in loop-diagram corrections to baryon number violating 
bosonic decays. Since it is necessary that the particles in the loop also undergo S-violating 
decays, the relevant particles are the X, Y, and H 3 bosons in the case of SU(5). In that 
case, CP violation is due to the complex phases of the Yukawa couplings h v and h D which 
cannot be reabsorbed by field redefinition. At tree-level these phases do not give any con- 
tribution to the baryon asymmetry and at the one-loop level the asymmetry is proportional 
to Im Tr (h\jhuh^ D hj^j = 0, where the trace is over generation indices [|19j. The problem of 
too little CP violation in SU(5) may be solved by further complicating the Higgs sector. 
For example, adding a Higgs in the 45 representation of SU(5) leads to an adequate baryon 
asymmetry for a wide range of the parameters ||48|| . In the case of SO(10), CP violation may 
be provided by the complex Yukawa couplings between the right-handed and the left-handed 
neutrinos and the scalar Higgs. 

2. CP violation in the CKM matrix of the Standard Model. 

Since only the left-handed electron is SU(2)l gauge coupled, C is maximally broken in 
the SM. Moreover, CP is known not to be an exact symmetry of the weak interactions. 
This is seen experimentally in the neutral Kaon system through K Q , K mixing. At present 
there is no accepted theoretical explanation of this. However, it is true that CP violation 
is a natural feature of the standard electroweak model. There exists a charged current 
which, in the weak interaction basis, may be written as Cw = (g / ' V / 2)^l7 m -DlH / a1 + h.c, 
where Ul = (u, c, t, . . .) L and D L = (d, s,b, . . Now, the quark mass matrices may be 
diagonalized by unitary matrices , V^, V®, V£ via 

diag(m u ,m c ,m t ,...) = Vg ' M U V% , (2.10) 
diag(m d , m S) m b) . . .) = V° M D V° . (2.11) 

Thus, in the basis of quark mass eigenstates, the charged current may be rewritten as 
C w = {g/V2)U' L KYD' L W^ + h.c, where U' L = V?U L and D' L = Vj°D L . The matrix 
K, defined by K = V^iV^y , is referred to as the Kobayashi-Maskawa (KM) quark mass 
mixing matrix. For three generations, as in the SM, there is precisely one independent 
nonzero phase 5 which signals CP- violation. While this is encouraging for baryogenesis, it 
turns out that this particular source of CP violation is not strong enough. The relevant 
effects are parameterized by a dimensionless constant which is no larger than 10 -20 . This 
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appears to be much too small to account for the observed BAU and, thus far, attempts 
to utilize this source of CP violation for electroweak baryogenesis have been unsuccessful. 
In light of this, it is usual to extend the SM in some minimal fashion that increases the 
amount of CP violation in the theory while not leading to results that conflict with current 
experimental data. 



3. CP violation in super symmetric models 



The most promising and well-motivated framework incorporating CP violation beyond 



the SM seems to be supersymmetry f4"5]4lj|. Let us consider the MSSM superpotential 

W = iiHiH 2 + h u H 2 Qu c + h d H l Qd c + h e H x Le c , (2.12) 

where we have omitted the generation indices. Here H\ and H 2 represent the two Higgs 
superfields, Q, u and d are the quark doublet and the up-quark and down-quark singlet 
superfields respectively and L and e c are the left-handed doublet and right-handed lepton 
singlet superfields. 

The lepton Yukawa matrix h e can be always taken real and diagonal while h u and h d con- 
tain the KM phase. There are four new sources for explicit CP violating phases, all arising 
from parameters that softly break supersymmetry. These are (i) the trilinear couplings 

T u H 2 Qu c + V d H 1 Qd c + Y e R x LZ c + h.c. , (2.13) 

where we have defined 

p(u,d,e) = J^u,d,e) . ^(w,d,e) ^ (2.14) 

and the tildes stand for scalar fields, (ii) the bilinear coupling in the Higgs sector nBH\H 2 , 
(iii) the gaugino masses M (one for each gauge group), and (iv) the soft scalar masses m. 
Two phases may be removed by redefining the phase of the superfield H 2 in such a way that 
the phase of \x is opposite to that of B. The product \iB is therefore real. It is also possible 
to remove the phase of the gaugino masses M by an R symmetry transformation. The latter 
leaves all the other supersymmetric couplings invariant, and only modifies the trilinear ones, 
which get multiplied by exp(— (J)m) where 0m is the phase of M. The phases which remain 
are therefore 



(f) A = arg(AM) and <^ = -arg(B). (2.15) 

These new phases <pA and M may be crucial for the generation of the baryon asymmetry. 

In the MSSM, however, there are other possible sources of CP violation. In fact, when 
supersymmetry breaking occurs, as we know it must, the interactions of the Higgs fields 
H\ and H 2 with charginos, neutralinos and stops (the superpartners of the charged, neutral 
gauge bosons, Higgs bosons and tops, respectively) at the one-loop level, lead to a CP 
violating contribution to the scalar potential of the form 

V CP = \x{H x H 2 f + X^H^H^ + A 3 |#2| 2 #i# 2 + h.c. (2.16) 
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These corrections may be quite large at finite temperature j|£],|5(J. One may write the 
Higgs fields in unitary gauge as Hi = 0) T and H 2 = (0, y? 2 e l6l ) T , where <pi, <p 2 , 9 are 
real, and 9 is the CP-odd phase. The coefficients Ai >2)3 determine whether the CP-odd 
field 9 is non zero, signalling the spontaneous violation of CP in a similar way to what 
happens in a generic two-Higgs model [^TJ. If, during the evolution of the early universe, a 
variation of the phase 9 is induced, this CP breakdown may bias the production of baryons 
through sphaleron processes in electroweak baryogenesis models. The attractive feature of 
this possibility is that, when the universe cools down, the finite temperature loop corrections 
to the A-couplings disappear, and the phase 9 relaxes to zero. This is turn implies that, 
unlike for scenarios utilizing <pA and one need not worry about present experimental 
constraints from the physics of CP violation. A recent study |52| has shown that when the 
constraints from experiment and the strength of the transition are taken into account, the 
relevant region of the parameter space for spontaneous CP violation in the MSSM at finite 
temperature is rather restricted. Nevertheless, in this small region, perturbative estimates 
need not be reliable, and non-perturbatively the region might be slightly larger (or smaller). 
More precise perturbative studies and/or lattice studies of the non-perturbative corrections 
are certainly needed. 

Finally, in scenarios in which the BAU is generated from the coherent production of a 
scalar condensate along a flat direction of the supersymmetric extension of the SM, CP- 
violation is present in the nonrenormalizable superpotentials which lift the flat directions at 
large field values. We shall have more to say about this is section ([V|). 



D. Departure from Thermal Equlibrium 

1. The out- of- equilibrium decay scenario 

Scenarios in which the third Sakharov condition is satisfied due to the presence of su- 
perheavy decaying particles in a rapidly expanding universe, generically fall under the name 
of out-of-equilibrium decay mechanisms |6|,[19[]. The underlying idea is fairly simple. If the 
decay rate Tx of the superheavy particles X at the time they become nonrelativistic [i.e. at 
the temperature T ~ Mx) is much smaller than the expansion rate of the universe, then the 
X particles cannot decay on the time scale of the expansion and so they remain as abundant 
as photons for T < Mx- In other words, at some temperature T > Mx, the superheavy 
particles X are so weakly interacting that they cannot catch up with the expansion of the 
universe and they decouple from the thermal bath while they are still relativistic, so that 
n x ~ n -y ~ T 3 at the time of decoupling. Therefore, at temperature T ~ M x , they popu- 
late the universe with an abundance which is much larger than the equilibrium one. This 
overbundance is precisely the departure from thermal equilibrium needed to produce a final 
nonvanishing baryon asymmetry when the heavy states X undergo B and CP violating de- 
cays. The out-of-equilibrium condition requires very heavy states: Mx ^ (10 15 — 10 16 ) GeV 
and M x *S (10 10 — 10 16 ) GeV, for gauge and scalar bosons, respectively [flQ[1 , if these heavy 
particles decay through renormalizable operators. 
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E. The Electroweak Phase Transition Scenario 



At temperatures around the electroweak scale, the expansion rate of the universe in ther- 
mal units is small compared to the rate of baryon number violating processes. This means 
that the equilibrium description of particle phenomena is extremely accurate at electroweak 
temperatures. Thus, baryogenesis cannot occur at such low scales without the aid of phase 
transitions and the question of the order of the electroweak phase transition becomes central. 

If the EWPT is second order or a continuous crossover, the associated departure from 
equilibrium is insufficient to lead to relevant baryon number production E3. This means 
that for EWBG to succeed, we either need the EWPT to be strongly first order or other 
methods of destroying thermal equilibrium to be present at the phase transition. 

Any thermodynamic quantity that undergoes such a discontinuous change at a phase 
transition is referred to as an order parameter, denoted by (p. For a first order transition 
there is an extremum at <p = which becomes separated from a second local minimum 
by an energy barrier. At the critical temperature T = T c both phases are equally favored 
energetically and at later times the minimum at tp 7^ becomes the global minimum of the 
theory. The dynamics of the phase tansition in this situation is crucial to most scenarios of 
electroweak baryogenesis. The essential picture is that around T c quantum tunneling occurs 
and nucleation of bubbles of the true vacuum in the sea of false begins. Initially these 
bubbles are not large enough for their volume energy to overcome the competing surface 
tension and they shrink and disappear. However, at a particular temperature below T c , 
bubbles just large enough to grow nucleate. These are termed critical bubbles, and they 
expand, eventually filling all of space and completing the transition. As the bubble walls 
pass each point in space, the order parameter changes rapidly, as do the other fields, and 
this leads to a significant departure from thermal equilibrium. Thus, if the phase transition 
is strongly enough first order it is possible to satisfy the third Sakharov criterion in this way. 

The precise evolution of critical bubbles in the electroweak phase transition is a crucial 
factor in determining which regimes of electroweak baryogenesis are both possible and effi- 
cient enough to produce the BAU. (See for possible obstacles to the standard picture) 
We discuss this, and the issue of the order of the phase transition in detail in section [TV]. 

For the bubble wall scenario to be successful, there is a further criterion to be satisfied. 
As the wall passes a point in space, the Higgs fields evolve rapidly and the Higgs VEV 
changes from (0) = in the unbroken phase to 



in the broken phase. Here, v(T) is the value of the order parameter at the symmetry breaking 
global minimum of the finite temperature effective potential. Now, CP violation and the 
departure from equilibrium occur while the Higgs field is changing. Afterwards, the point is 
in the true vacuum, baryogenesis has ended, and baryon number violation is exponentially 
supressed. Since baryogenesis is now over, it is imperative that baryon number violation be 
negligible at this temperature in the broken phase, otherwise any baryonic excess generated 
will be equilibrated to zero. Such an effect is known as washout of the asymmetry and the 
criterion for this not to happen may be written as 




(2.17) 




c 



(2.18) 
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Although there are a number of nontrivial steps that lead to this simple criterion, (|2.18|) is 
traditionally used to ensure that the baryon asymmetry survives after the wall has passed. 
It is necessary that this criterion be satisfied for any electroweak baryogenesis scenario to 
be successful. 



1. Defect-Mediated Baryogenesis 



A natural way to depart from equilibrium is provided by the dynamics of topological 
defects. Topological defects are regions of trapped energy density which can remain after 
a cosmological phase transition if the topology of the vacuum of the theory is nontrivial. 
Typically, cosmological phase transitions occur when a symmetry of a particle physics the- 
ory is spontaneously broken. In that case, the cores of the topological defects formed are 
regions in which the symmetry of the unbroken theory is restored. If, for example, cosmic 
strings are produced at the GUT phase transition, then the decays of loops of string act 
as an additional source of superheavy bosons which undergo baryon number violating de- 
cays. When defects are produced at the TeV scale, a detailed analysis of the dynamics of 
a network of these objects shows that a significant baryon to entropy ratio can be gener- 
ated [^T|,0,^,^,^,^] if the electroweak symmetry is restored around such a higher scale 
ordinary defect [p^-|69|. Although a recent analysis has shown that B- violation is highly 
inefficient along nonsuper conducting strings ||70[| , there remain viable scenarios involving 
other ordinary defects, superconducting strings or defects carrying baryon number [|7T],[72|. 



III. RECENT DEVELOPMENTS IN GUT BARYOGENESIS 



As we have seen, in GUTs baryon number violation is natural since quarks and leptons 
lie in the same irreducible representation of the gauge group. The implementation of the 
out-of-equilibrium scenario in this context goes under the name of GUT baryogenesis. The 
basic assumption is that the superheavy bosons were as abundant as photons at very high 
temperatures T J> Mj. This assumption is questionable if the heavy X particles are the 
gauge or Higgs bosons of Grand Unification, because they might have never been in thermal 
equilibrium in the early universe. Indeed, the temperature of the universe might have been 
always smaller than Mqut and, correspondingly, the thermally produced X bosons might 
never have been as abundant as photons, making their role in baryogenesis negligible. All 
these considerations depend crucially upon the thermal history of the universe and deserve 
a closer look. 



A. Inflation, reheating and GUT baryogenesis: the old wisdom 



The flatness and the horizon problems of the standard big bang cosmology are elegantly 
solved if, during the evolution of the early universe, the energy density happened to be dom- 
inated by some form of vacuum energy, causing comoving scales grow quasi-exponentially 
76^791 . The vacuum energy driving this inflation is generally assumed to be associated 



with the potential V{(j>) of some scalar field 0, the inflaton, which is initially displaced from 
the minimum of its potential. As a by-product, quantum fluctuations of the inflaton field 
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may be the seeds for the generation of structure and the fluctuations observed in the cosmic 



microwave background radiation ||79|| . Inflation ended when the potential energy associated 
with the inflaton field became smaller than the kinetic energy of the field. By that time, 
any pre-inflationary entropy in the universe had been inflated away, and the energy of the 
universe was entirely in the form of coherent oscillations of the inflaton condensate around 
the minimum of its potential. Now, we know that somehow this low-entropy cold universe 
must be transformed into a high-entropy hot universe dominated by radiation. The process 
by which the energy of the inflaton field is transferred from the inflaton field to radiation has 
been dubbed reheating. The simplest way to envision this process is if the comoving energy 
density in the zero mode of the inflaton decays perturbatively into normal particles, which 
then scatter and thermalize to form a thermal background pOlJSll . It is usually assumed 
that the decay width of this process is the same as the decay width of a free inflaton field. 
Of particular interest is a quantity known as the reheat temperature, denoted as T r . This 
is calculated by assuming an instantaneous conversion of the energy density in the inflaton 
field into radiation when the decay width of the inflaton energy, T^, is equal to H, the 
expansion rate of the universe. This yields 



T r = sjT^rrip. (3.1) 

There are very good reasons to suspect that GUT baryogenesis does not occur if this is the 
way reheating happens. First of all, the density and temperature fluctuations observed in 
the present universe require the inflaton potential to be extremely flat - that is -C M^. 
This means that the couplings of the inflaton field to the other degrees of freedom cannot be 
too large, since large couplings would induce large loop corrections to the inflaton potential, 
spoiling its flatness. As a result, T r is expected to be smaller than 10 14 GeV by several 
orders of magnitude. On the other hand, the unification scale is generally assumed to be 
around 10 16 GeV, and S-violating bosons should have masses comparable to this scale. 
Furthermore, even the light 5-violating Higgs bosons are expected to have masses larger 
than the inflaton mass, and thus it would be kinematically impossible to create them directly 
in decays, (p — > XX. 

One might think that the heavy bosons could be created during the stage of thermal- 
ization. Indeed, it has recently been shown that particles of mass ~ 10 4 bigger than the 
reheating temperature T r may be created by the thermalized decay products of the inflaton 
H82HS5|. However, there is one more problem associated with GUT baryogenesis in the old 
theory of reheating, namely the problem of relic gravitinos [^BJ. If one has to invoke su- 



persymmetry to preserve the flatness of the inflaton potential, it is mandatory to consider 
the cosmological implications of the gravitino - a spin- (3/2) particle which appears in the 
extension of global supersymmetry to supergravity. The gravitino is the fermionic super- 
partner of the graviton and has interaction strength with the observable sector - that is 
the standard model particles and their superpartners - inversely proportional to the Planck 
mass. The slow gravitino decay rate leads to a cosmological problem because the decay 
products of the gravitino destroy 4 He and D nuclei by photodissociation, and thus ruin the 
successful predictions of nucleosynthesis. The requirement that not too many gravitinos 
are produced after inflation provides a stringent constraint on the reheating temperature, 
T r <, (10 10 — 10 11 ) GeV |p6| . Therefore, if T r ~ M G ut, gravitinos would be abundant during 
nucleosynthesis and destroy the agreement of the theory with observations. However, if the 
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initial state after inflation was free from gravitinos, the reheating temperature is then too low 
to create superheavy X bosons that eventually decay and produce the baryon asymmetry. 



B. Inflation, reheating and GUT baryogenesis: the new wisdom 

The outlook for GUT baryogenesis has brightened recently with the realization that 
reheating may differ significantly from the simple picture described above fl87-9lH. In the first 



stage of reheating, called preheating [[37|| , nonlinear quantum effects may lead to extremely 
effective dissipative dynamics and explosive particle production, even when single particle 
decay is kinematically forbidden. In this picture, particles can be produced in a regime 
of broad parametric resonance, and it is possible that a significant fraction of the energy 
stored in the form of coherent inflaton oscillations at the end of inflation is released after 
only a dozen or so oscillation periods of the inflaton. What is most relevant for us is that 
preheating may play an extremely important role in baryogenesis and, in particular, for the 
Grand Unified generation of a baryonic excess. Indeed, it was shown in [0,0 that the 
baryon asymmetry can be produced efficiently just after the preheating era, thus solving 
many of the problems that GUT baryogenesis had to face in the old picture of reheating. 

The presence of a preheating stage is due to the fact that, for some parameter ranges, 
there is a new non-perturbative decay channel. Coherent oscillations of the inflaton field 
induce stimulated emissions of bosonic particles into energy bands with large occupancy 
numbers ||37|| . The modes in these bands can be understood as Bose condensates, and they 
behave like classical waves. 

The idea of preheating is relatively simple, the oscillations of the inflaton field induce a 
time-dependent mass for fields coupled to the inflaton. This leads to a mixing of positive 
and negative frequencies in the quantum state of these fields. For the sake of simplicity, let 
us focus on the case of chaotic inflation, with a massive inflaton <f), with quadratic potential 
V(4>) = |M?0 2 , ~ 10 13 GeV, coupled to a massive scalar field X via the quartic coupling 
g 2 (p 2 \X\ 2 . Let us also neglect the expansion of the universe. 

The evolution equation for the Fourier modes of the X field with momentum k is 

X k + u 2 k X k = 0, (3.2) 

with io\ = k 2 + M\ + g 2 <j) 2 (t). This Klein-Gordon equation may be cast in the form of a 
Mathieu equation 

X' k ' + [A(k) - 2q cos 2z]X k = 0, (3.3) 



where z = M^t and 



k 2 + M 2 x 



$2 



f^jjjj, (3-4) 

where $ is the amplitude and is the frequency of inflaton oscillations, defined through 
<f)(t) = $(t) sin(M^t). Notice that, at least initially when $ = cm p <> m p 
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^TTH ~ 9 2 c 2 ^ ~g 2 c 2 x 10 12 » 1 (3.5) 

and the resonance is broad. A crucial observation for baryogenesis is that particles with 
mass larger than that of the inflaton may be produced during preheating P2"| . Particle 
production occurs above the line A = 2q, and the width of the instability strip scales as q 1 ^ 2 
for large q, independent of Mx- The condition for broad resonance |87| , |92| A — 2q <, q 1 ^ 2 
becomes (k 2 + M 2 x )/Ml <, g{$/M+), which yields E\ = k 2 + M| <> g$Mf for the typical 
energy of X bosons produced in preheating. By the time the resonance develops to the full 
strength, $ 2 ~ 10~ 5 m 2 , and so the resulting estimate for the typical energy of particles at 
the end of the broad resonance regime for Ma, ~ 10 13 GeV is 



E x ~ ~ g 1/2 !0 15 GeV. (3.6) 

Supermassive X bosons can be produced by the broad parametric resonance for Ex > Mx, 
which leads to the condition M x < g 1/2 10 15 GeV |[7||)2f. Thus, for g 2 ~ 1 one would have 
copious production of X particles as heavy as 10 15 GeV, i.e., 100 times greater than the 
inflaton mass. 

The corresponding number density n k grows exponentially with time because these par- 
ticles belong to an instability band of the Mathieu equation (for a recent comprehensive 



review on preheating after chaotic inflation pi and references therein) 



X k oc e^ kM ^ ^n k (x e 2 ^ A V (3.7) 

where the parameter fi k depends upon the instability band and, in the broad resonance case, 
q > 1, it is ~ 0.2. 

These instabilities can be interpreted as coherent "particle" production with large occu- 
pancy numbers. One way of understanding this phenomenon is to consider the energy of 
these modes as that of a harmonic oscillator, 

E k = \X k \ 2 /2 + uo 2 k \X k \ 2 /2 = uo k {n k + 1/2) . (3.8) 

The occupancy number of level k can grow exponentially fast and these modes soon behave 
like classical waves. The parameter q during preheating determines the strength of the 
resonance. 

This is a significant departure from the old constraints of reheating. Production of X 
bosons in the old reheating picture was kinematically forbidden if < Mx, while in the 
new scenario it is possible because of coherent effects. Not only are particles produced out-of- 
equilibrium, thus satisfying one of the basic requirements to produce the baryon asymmetry, 
but also gravitinos may be shown to be produced in sufficiently small amounts 0| to avoid 
the gravitino problem we discussed earlier. 

Back-reaction, or scattering of X fluctuations off the zero mode of the inflaton field, 
limits the maximum magnitude of X fluctuations (X 2 ) P0 |. This in turn restricts the corre- 



sponding number density of created V-particles and therefore weakens the whole mechanism. 
Parametric resonance is also rendered less efficient when the X particles have a decay width 



T x , a crucial parameter in GUT baryogenesis. Using the methods developed in Refs. pq-90 



the problem of the production of massive, unstable X particles in the process of the inflaton 
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decay has been recently addressed |53| for the case of a quadratic inflaton potential. It was 
shown that the observed baryon asymmetry may be explained by GUT baryogenesis after 
preheating if 



where T = Tx/M^, and e is a dimensionless measure of the CP violation in the system. 
This may happen if |93j Fx ^ 10 -3 Mj, which is a relatively mild constraint. 

Grand Unified baryogenesis after preheating solves many of the serious drawbacks of 
GUT baryogenesis in the old theory of reheating where the production of superheavy states 
after inflation was kinematically forbidden. Moreover, the out-of-equilibrium condition is 
naturally achieved in this scenario, since the distribution function of the X-quanta generated 
at the resonance is far from a thermal distribution. This situation is considerably different 
from the one present in the thermal GUT scenario where superheavy particles usually de- 
couple from the thermal bath when still relativistic and then decay producing the baryon 
asymmetry. It is quite intriguing that out of all possible ways the parametric resonance may 
develop, Nature might have chosen only those ways without instantaneous thermalization 
and also with a successful baryogenesis scenario. 



C. Baryogenesis via leptogenesis 

Since the linear combination B — L is left unchanged by sphaleron transitions, the baryon 
asymmetry may be generated from a lepton asymmetry ||23|| . Indeed, sphaleron transition 
will reprocess any lepton asymmetry and convert (a fraction of) it into baryon number. This 
is because B + L must be vanishing and the final baryon asymmetry results to be B ~ —L. 
Once the lepton number is produced, thermal scatterings redistribute the charges. In the 
high temperature phase of the SM, the asymmetries of baryon number B and of B — L are 
therefore proportional: |I3 



/ 8nf + Anu \ , , . 

B=[ f — — )(B-L), 3.10 

\22nf + 13n H J 

where rifj is the number of Higgs doublets. In the SM as well as in its unified extension 
based on the group SU(5), B — L is conserved and no asymmetry in B — L can be gener- 
ated. However, adding right-handed Majorana neutrinos to the SM breaks B — L and the 
primordial lepton asymmetry may be generated by the out-of-equilibrium decay of heavy 
right-handed Majorana neutrinos iV£ (in the supersymmetric version, heavy scalar neutrino 
decays are also relevant for leptogenesis). This simple extension of the SM can be embedded 
into GUTs with gauge groups containing SO( 10). Heavy right-handed Majorana neutrinos 
can also explain the smallness of the light neutrino masses via the see-saw mechanism p5 



The relevant couplings are those between the right-handed neutrino, the Higgs doublet 
$ and the lepton doublet £l 

C = ~ L $KNl + hr L MN c L + h. C . (3.11) 
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The vacuum expectation value of the Higgs field ($) generates neutrino Dirac masses ran = 
h v (<&), which are assumed to be much smaller than the Majorana masses M. When the 
Majorana right-handed neutrinos decay into leptons and Higgs scalars, they violate the 
lepton number since right-handed neutrino fermionic lines do not have any preferred arrow 

N C L -> $ + £, 

N C L ^ $ + £ (3.12) 

The interference between the tree-level decay amplitude and the absorptive part of the one- 
loop vertex leads to a lepton asymmetry of the right order of magnitude to explain the 
observed baryon asymmetry and has been discussed extensively in the literature [p3|,|96|-98 



Recently, much attention has been paid to the effects of finite temperature on the CP 
violation ||99|| , and the contribution to the lepton asymmetry generated by the tree- level 



graph and the absorptive part of the one-loop self-energy. In particular, it has been observed 
that CP violation may be considerably enhanced if two heavy right-handed neutrinos are 
nearly degenerate in mass [ 100|| . It is also important to remember that large lepton number 



violation at intermediate temperatures may potentially dissipate away the baryon number 
in combination with the sphaleron transitions. Indeed, AL = 2 interactions of the form 

777 

"" ^ L M + h.c, (3.13) 



where m u is the mass of the light left-handed neutrino, are generated through the exchange 
of heavy right-handed neutrinos. The rate of lepton number violation induced by this 
interaction is therefore Yl ~ (m 2 /($) 4 )T 3 . The requirement of harmless letpon number 
violation, Yl < H imposes an interesting bound on the neutrino mass 

m ^ 4eV (lo^f /2 - (314) 

where Tx = Min {Tb~l, 10 12 GeV}, Tb-l is the temperature at which the B — L number 
production takes place, and ~ 10 12 GeV is the temperature at which sphaleron transitions 
enter in equilibrium. One can also reverse the argument and study leptogenesis assuming 
a similar pattern of mixings and masses for leptons and quarks, as suggested by SO (10) 
unification [ 101 1 . This implies that B — L is broken at the unification scale ~ 10 16 GeV, 



if m„ ~ 3 x 10 3 eV as preferred by the MSW explanation of the solar neutrino deficit 
1021,110311. 



IV. ELECTROWEAK BARYOGENESIS 

Scenarios in which anomalous baryon number violation in the electroweak theory im- 
plements the first Sakharov criterion are refered to as electroweak baryogenesis. Typically 
these scenarios require a strongly first order phase transition, either in the minimal standard 
model, or in a modest extension. Another possibility is that the out of equilibrium require- 
ment is achieved through TeV scale topological defects, also present in many extensions of 
the electroweak theory. Finally, a source of CP-violation beyond that in the CKM matrix 
is usually involed. 
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Historically, the ways in which baryons may be produced as a bubble wall, or phase 
boundary, sweeps through space, have been separated into two categories. 



1. local baryogenesis: baryons are produced when the baryon number violating processes 
and CP violating processes occur together near the bubble walls. 

2. nonlocal baryogenesis: particles undergo CP violating interactions with the bubble 
wall and carry an asymmetry in a quantum number other than baryon number into 
the unbroken phase region away from the wall. Baryons are then produced as baryon 
number violating processes convert the existing asymmetry into one in baryon number. 

In general, both local and nonlocal baryogenesis will occur and the BAU will be the sum 
of that generated by the two processes. However, if the speed of the wall is greater than the 
sound speed in the plasma, then local baryogenesis dominates [|104j| . In other cases, nonlocal 
baryogenesis is usually more efficient and we shall focus on that here. 

Nonlocal baryogenesis typically involves the interaction of the bubble wall with the var- 
ious fermionic species in the unbroken phase. The main picture is that as a result of CP 
violation in the bubble wall, particles with opposite chirality interact differently with the 
wall, resulting in a net injected chiral flux. This flux thermalizes and diffuses into the 
unbroken phase where it is converted to baryons. In this section, for definiteness when de- 
scribing these effects, we assume that the CP violation arises because of a two-Higgs doublet 
structure. 

The chiral asymmetry which is converted to an asymmetry in baryon number is carried 
by both quarks and leptons. However, the Yukawa couplings of the top quark and the 
r-lepton are larger than those of the other quarks and leptons respectively. Therefore, it 
is reasonable to expect that the main contribution to the injected asymmetry comes from 
these particles and to neglect the effects of the other particles. 

When considering nonlocal baryogenesis it is convenient to write the equation for the 
rate of production of baryons in the form [[10 5| 



dn B n f T sp {T) x . 

dt ~ 2T 2 -^' 1 ' 

where the rate per unit volume for electroweak sphaleron transitions is given by ( |2.9| ). Here, 
nf is again the number of families and fii is the chemical potential for left handed particles 
of species i. The crucial question in applying this equation is an accurate evaluation of the 
chemical potentials that bias baryon number production. 

There are typically two distinct calculational regimes that are appropriate for the treat- 
ment of nonlocal effects in electroweak baryogenesis. Which regime is appropriate depends 
on the particular fermionic species under consideration. 



1. The thin wall regime 



The thin wall regime [11061 , 11071 , 110511 applies if the mean free path I of the fermions being 
considered is much greater than the thickness L w of the wall, i.e. if 
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In this case we may neglect scattering effects and treat the fermions as free particles in their 
interactions with the wall. 

In this regime, in the rest frame of the bubble wall, particles see a sharp potential barrier 
and undergo CP violating interactions with the wall due to the gradient in the CP odd 
Higgs phase. As a consequence of CP violation, there will be asymmetric reflection and 
transmission of particles, thus generating an injected current into the unbroken phase in 
front of the bubble wall. As a consequence of this injected current, asymmetries in certain 
quantum numbers will diffuse both behind and in front of the wall due to particle interactions 
and decays ||106| , |107| , |105|| . In particular, the asymmetric reflection and transmission of left 
and right handed particles will lead to a net injected chiral flux from the wall. However, 
there is a qualitative difference between the diffusion occurring in the interior and exterior 
of the bubble. 

Exterior to the bubble the electroweak symmetry is restored and weak sphaleron tran- 
sitions are unsuppressed. This means that the chiral asymmetry carried into this region by 
transport of the injected particles may be converted to an asymmetry in baryon number by 
sphaleron effects. In contrast, particles injected into the phase of broken symmetry interior 
to the bubble may diffuse only by baryon number conserving decays since the electroweak 
sphaleron rate is exponentially suppressed in this region. Hence, we concentrate only on 
those particles injected into the unbroken phase. 

The net baryon to entropy ratio which results via nonlocal baryogenesis in the case of thin 
walls has been calculated in several different analyses ||106| , p.07|| and | 105 |. In the following 



we give a brief outline of the logic of the calculation, following ||105||. The baryon density 



produced is given by ( |4. 1|) in terms of the chemical potentials for left handed particles. 
These chemical potentials are a consequence of the asymmetric reflection and transmission 
off the walls and the resulting chiral particle asymmetry. Baryon number violation is driven 
by the chemical potentials for left handed leptons or quarks. To be concrete, we focus on 
leptons [|105|| (for quarks see e.g. [|106|| ). If there is local thermal equilibrium in front of the 
bubble walls then the chemical potentials fii of particle species i are related to their number 
densities rii by 

T 2 

rii = j^ki^i , (4.3) 

where hi is a statistical factor which equals 1 for fermions and 2 for bosons. In deriving this 
expression, it is important |108| to correctly impose the constraints on quantities which are 



conserved in the region in front of and on the wall. 

Using the above considerations, the chemical potential \xl for left handed leptons can 
be related to the left handed lepton number densities L^. These are in turn determined by 
particle transport. The source term in the diffusion equation is the flux J resulting from 
the asymmetric reflection and transmission of left and right handed leptons off the bubble 
wall. 

For simplicity assume a planar wall. If \p z \ is the momentum of the lepton perpendic- 
ular to the wall (in the wall frame), the analytic approximation used in ||105|| allows the 



asymmetric reflection coefficients for lepton scattering to be calculated in the range 
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m t < \p z \ < m h ~ — , (4.4) 
where mi and are the lepton and Higgs masses, respectively, and results in 

2 

K L ^ R - K R ^ L ~ 2M CP -^ . (4.5) 

m h \p z \ 

The corresponding flux of left handed leptons is 

v w mfm h A6 C p 



4vr 2 



(4.6) 



where v w is the velocity of the bubble wall. Note that in order for the momentum interval 
in ( |4.4j) to be non-vanishing, the condition miL w < 1 needs to be satisfied. 

The injected current from a bubble wall will lead to a "diffusion tail" of particles in front 
of the moving wall. In the approximation in which the persistence length of the injected 
current is much larger than the wall thickness we may to a good approximation model it as 
a delta function source and search for a steady state solution. In addition, assume that the 
decay time of leptons is much longer than the time it takes for a wall to pass so that we may 
neglect decays. Then the diffusion equation for a single particle species is easily solved by 

- j °- e : ; i . (4.7) 



with the diffusion root 



A^ = ^. M.N) 



L 

where Dl is the diffusion constant for leptons, and £ L that is called the persistence length 
of the current in front of the bubble wall. 

Note that in this approximation the injected current does not generate any perturbation 
behind the wall. This is true provided £ L 3> L w is satisfied. If this inequality is not true, 
the problem becomes significantly more complex [|105|| . 

In the massless approximation the chemical potential fi^ can be related to by 

t*L = ^L L (4.9) 

(for details see [ |105|| ). Inserting the sphaleron rate and the above results for the chemical 
potential /i into ([4.1|), and using 1/Dl — 8a^/T, the final baryon to entropy ratio becomes 

^O.2a. ; „ (j TVA0 o A(^) 2 ^f ■ (4.10) 
s v w \ 1 J 1 D L 

where we have assumed that sphalerons do not equilibrate in the diffusion tail. 

Now consider the effects of top quarks scattering off the advancing wall | |106| , |107| |. Several 
effects tend to decrease the contribution of the top quarks relative to that of tau leptons. 
Firstly, for typical wall thicknesses the thin wall approximation does not hold for top quarks. 
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This is because top quarks are much more strongly interacting than leptons and so have a 
much shorter mean free path. An important effect is that the diffusion tail is cut off in front 



of the wall by strong sphalerons ||109| , |110|| . There is an anomaly in the quark axial vector 
current in QCD. This leads to chirality non- conserving processes at high temperatures. 
These processes are relevant for nonlocal baryogenesis since it is the chirality of the injected 
current that is important in that scenario. In an analogous expression to that for weak 
sphalerons, we may write the rate per unit volume of chirality violating processes due to 
strong sphalerons in the unbroken phase as 

T s = K s {a s Tf , (4.11) 



where k s is a dimensionless constant [ 111 ] . Note that the uncertainties in k s are precisely the 
same as those in k defined in (|2.9| ). As such, k s could easily be proportional to a s , in analogy 
with ( |2.9| ), perhaps with a logarithmic correction. These chirality-changing processes damp 
the effect of the injected chiral flux and effectively cut off the diffusion tail in front of the 
advancing bubble wall. Second, the diffusion length for top quarks is intrinsically smaller 
than that for tau leptons, thus reducing the volume in which baryogenesis takes place. 
Although there are also enhancement factors, e.g. the ratio of the squares of the masses 
m^/rn^, it seems that leptons provide the dominant contribution to nonlocal baryogenesis. 



2. The thick wall regime 

Now we turn to the thick wall, or adiabatic, regime. This is relevant if the mean free path 
of the particles being considered is smaller than the width of the wall, i ^ L w . When the 
walls are thick, most interactions within the wall will be almost in thermal equilibrium. The 
equilibrium is not exact because some interactions, in particular baryon number violation, 
take place on a time scale much slower than the rate of passage of the bubble wall. These 
slowly-varying quantities are best treated by the method of chemical potentials. The basic 
idea is to produce asymmetries in some local charges which are (approximately) conserved by 
the interactions inside the bubble walls, where local departure from thermal equilibrium is 
attained. These local charges will then diffuse into the unbroken phase where baryon number 
violation is active thanks to the unsuppressed sphaleron transitions. The latter convert the 
asymmetries into a baryon asymmetry. Therefore, one has to i) identify those charges 
which are approximately conserved in the symmetric phase, so that they can efficiently 
diffuse in front of the bubble where baryon number violation is fast, ii) compute the CP 
violating currents of the plasma induced inside the bubble wall and Hi) solve a set of coupled 
differential diffusion equations for the local particle densities, including the CP violating 
sources. 

There are a number of possible CP violating sources. To be definite, consider the example 



where CP violation is due to a CP odd phase 9 in the two-Higgs doublet model ||112|| . This 
goes generically under the name of spontaneous baryogenesis. In order to explicitly see 
how 9 couples to the fermionic sector of the theory (to produce baryons) we may remove 
the ^-dependence of the Yukawa couplings arising from the Higgs terms. We do this by 



performing an anomaly- free hypercharge rotation on the fermions [112 1 , inducing a term 



in the Lagrangian density of the form Ccp oc d^Jy, where J Y is the fermionic part of 



21 



the hypercharge current. Therefore, a nonvanishing 9 provides a preferential direction for 
the production of quarks and leptons; in essence a chemical potential fiB oc 9 (in the rest 
frame of the thermal bath) for baryon number. Of course, strictly speaking, this is not a 
chemical potential since it arises dynamically rather than to impose a constraint. For this 
reason, the quantity fiB is sometimes referred to as a "charge potential". A more complete 
field theoretic approach to the computation of particle currents on a space-time dependent 
and CP violating Higgs background was provided in [ |113|J114|| where it was shown that 
fermionic currents arise at one loop, and are proportional to (hv(T c )/7rT) 2 9, where h is 
the appropriate Yukawa coupling. The relevant baryon number produced by spontaneous 
baryogenesis is then calculated by integrating 
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Initially, spontaneous baryogenesis was considered as an example of local baryogenesis. How- 
ever, it has become clear that diffusion effects can lead to an appreciable enhancement of 
the baryon asymmetry produced by this mechanism [ 115 - 11^] ■ 

An alternative way of generating a source for the diffusion equation was suggested by 



Joyce et al. ||118|| . The essential idea is that there exists a purely classical chiral force that 
acts on particles when there is a CP violating field on the bubble wall. The Lagrangian for 
a fermion ip in the background of a bubble wall in the presence of the CP-odd phase 9 is 



c = a + 



2 v\ + v: 



^d^9 ip - mi/jip , 



(4.13) 



where V\ and v-i are the VEVs of the Higgs fields and m is the fermion mass. As usual we 
assume a planar wall propagating at speed v in the ^-direction. Further, if interactions of 
the Higgs and gauge fields with the wall have reached a stationary state, then these fields are 
functions only of z — vt. Then, in the wall rest frame plane- wave fermion solutions ip oc e~ ip ' x 
have the dispersion relation 
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(4.14) 



where ± corresponds to s z = ±1/2, with s z the z-component of the spin. 

In the WKB approximation, the Hamilton equations derived from this relation allow 
calculation of the acceleration of a WKB wave-packet (i.e. the chiral force). This force acts 
as a potential well that induces an excess of chiral charge on the wall. In the diffusion tail in 
front of the wall, there exists a corresponding deficit of chiral charge. This acts as a chemical 
potential which sources the production of baryon number. 

This calculation is performed by following a fluid approach to solving the associated 
Boltzmann equation, By concentrating on particles with \p z \ ~ T ^> m, the particle and 
antiparticle excitations ( |4.14| ) may be treated as separate fluids in the WKB approximation. 
This permits an analytic solution to the diffusion equation which yields a relation between 
the chemical potential on the wall for the fermion species considered and the chiral force. 
Integrating this in front of the wall yields a chemical potential that may then be fed into 
(|1.1|) to give an approximate answer for the BAU. 



22 



To summarize, the effects of transport in the plasma external to the expanding bubble 
allow baryon violating transitions in the unbroken phase to transform charge asymmetries 
produced on the wall into a baryon asymmetry. This means that in the case of nonlocal 
baryogenesis we do not need to rely on baryon number violating processes occurring in 
the region where the Higgs fields are changing. The diffusion equation approach to the 
problem of nonlocal baryogenesis has been very successful. In the thin wall case, it is a valid 
approximation to assume that the source for the diffusion equation is essentially a 5-function 
on the wall. This is because one may ignore the effects of particle scattering in this picture. 
However, in the case of thick walls, significant particle scattering occurs and as a result it is 
necessary to consider sources for the diffusion equations that extend over the wall. 



A. Electroweak baryogenesis in the MSSM 

As we mentioned, the most promising and well- motivated framework incorporating CP 
violation beyond the SM and an enhanced electroweak phase transition seems to be super- 
symmetry. Electroweak baryogenesis in the framework of the MSSM has attracted much 
attention in the past years, with particular emphasis on the strength of the phase transi- 
tion [TTS-O^] and the mechanism of baryon number generation ||1 23LII 24j,pT|,|l 2 3- 1 28 . 



1. The strength of the electroweak phase transition in the MSSM 



The behavior of the electroweak phase transition in the minimal supersymmetric stan- 
dard model is dependent on the mass of the lightest Higgs particle, and the mass of the 
top squark. Recent analytical [p.29| - |137| and lattice computations [|138| , |140|J139|J141|| have 
revealed that the phase transition can be sufficiently strongly first order in the presence of a 
top squark lighter than the top quark. In order to naturally suppress contributions to the p- 
parameter, and hence preserve a good agreement with precision electroweak measurements 
at LEP, the top squark should be mainly right handed. This can be achieved if the left 
handed stop soft supersymmetry breaking mass tuq is much larger than Mz- For moderate 
mixing, the lightest stop mass is then approximately given by 
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(4.15) 



where A t = A t — fi*/ tan (3 is the particular combination appearing in the off-diagonal terms 
of the left-right stop squared mass matrix and my is the soft supersymmetry breaking mass 
parameter of the right handed stop. 

The preservation of the baryon number asymmetry requires the order parameter 
(4>(T C ))/T C to be larger than one, see Eq. ( |2.18|) . This quantity is bounded from above 
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where m% = mt(mt) is the on-shell running top quark mass in the MS scheme, rrih is the 
lightest Higgs boson mass and v = Jvf + v\ is the vev of the Higgs fields. The first term 



on the right hand side of expression (4.16) is the Standard Model contribution 
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(4.17) 



and the second term is the contribution that would be obtained if the right handed stop 
plasma mass vanished at the critical temperature (see Eq. ( |4.18| )). The difference between 
the SM and the MSSM is that light stops may give a large contributions to the effective 
potential in the MSSM and therefore overcome the Standard Model constraints. The stop 
contribution strongly depends on the value of mfj, which must be small in magnitude, and 
negative, in order to induce a sufficiently strong first order phase transition. Indeed, large 
stop contributions are always associated with small values of the right handed stop plasma 



mass 
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(4.19) 



is the finite temperature self-energy contribution to the right-handed squarks. Moreover, 
the trilinear mass term, A t , must satisfy Af <C mg in order to avoid the suppression of the 
stop contribution to (<f>(T c ))/T c . Note that, although large values of mu, of order of the 
critical temperature, are useful to get a strongly first order phase transition, mu is bounded 



from above to avoid the appearance of dangerous charge and color breaking minima [ 136|| . 

As is clear from ([4.16 ), in order to obtain values of (0(T c ))/T c larger than one, the Higgs 
mass must take small values, close to the present experimental bound. For small mixing 
and large m^, the one-loop Higgs mass has a very simple form 
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where m~ ~ m 2 q + mf. Hence, small values of tan/3 are preferred. The larger the left 
handed stop mass, the closer to one tan/3 must be. This implies that the left handed stop 
effects are likely to decouple at the critical temperature, and hence that mq mainly affects 
the baryon asymmetry through the resulting Higgs mass. A detailed analysis, including all 
dominant two-loop finite temperature corrections to the Higgs effective potential and the 
non-trivial effects arising from mixing in the stop sector, has been recently performed ||136 



This analysis has identified the region of parameter space for which MSSM electroweak 
baryogenesis can happen. Taking into account the experimental bounds as well as the 
requirement of avoiding dangerous color breaking minima, it was found that the lightest 
Higgs should be lighter than about 105 GeV, while the stop mass may be close to the 
present experimental bound and must be smaller than, or of order of, the top quark mass. 
This lower bound has been essentially confirmed by lattice simulations | 141|| , providing a 
motivation for the search for Higgs and stop particles at the LEP and Tevatron colliders. 



24 



2. Baryon number generation in the MSSM 



As we have seen, the MSSM contains additional sources of CP-violation besides the 
CKM matrix phase. These new phases are essential for the generation of the baryon number 
since large CP violating sources may be locally induced by the passage of the bubble wall 
separating the broken from the unbroken phase during the electroweak phase transition. 
The new phases appear in the soft supersymmetry breaking parameters associated with the 
stop mixing angle and the gaugino and neutralino mass matrices. However, large values of 
the stop mixing angle are strongly restricted in order to preserve a sufficiently strong first 
order electroweak phase transition. Therefore, an acceptable baryon asymmetry may only be 
generated from the stop sector through a delicate balance between the values of the different 
soft supersymmetry breaking parameters contributing to the stop mixing parameter, and 
their associated CP violating phases [124]. As a result, the contribution to the final baryon 



asymmetry from the stop sector turns out to be negligible. On the other hand, charginos and 
neutralinos may be responsible for the observed baryon asymmetry if a large CP violating 
source for the higgsino number 7^ is induced in the bubble wall. This can happen for 



sufficiently large values of the phase of the parameter [i [|124| , |128| 

The first step in the computation of the baryon asymmetry is to solve a set of coupled 
differential equations describing the effects of diffusion, particle number changing reactions, 
and CP violating source terms ||123| . |124| , [L28|| . Major simplifications of the diffusion equations 
take place when neglecting all the couplings except for gauge interactions and the top Yukawa 
coupling. Neglecting the weak sphalerons (in the first step) allows us to ignore leptons in the 
diffusion equations, and turns out to be a good approximation when computing Higgs and 
quark densities. A further major simplification is to assume that the top Yukawa and the 
strong sphaleron interactions are fast, thus reducing the diffusion equations to an equation 
governing the baryon asymmetry tib- We write the integral of the CP violating source as 

D A+ Jo 

where 



1 t°° 

_4 K =— / du~,~ H (u)e-^" , (4.21) 



+ 2D y ' 



and T is the effective decay constant governing the particle number changing reaction ||123 
Then, the solution for the baryon to entropy ratio is 

2- = - 9 (fc)^P , (4.23) 
s vis 

where g{kj) is a numerical coefficient depending upon the light degrees of freedom present 
in the thermal bath, and D is an effective diffusion constant. CP violating densities are 
therefore non zero for a time t ~ D/v^ and the assumptions leading to the above expressions 
are valid as long as top Yukawa and the strong sphaleron interactions are faster than ~ D/v^ 
J128| . |124| |. On larger time scales, diffusion takes over, washing out the sources. 

The next step is the evaluation of the CP violating source for higgsino number. In this 
case, non-equilibrium quantum field theory provides the necessary tools to write down a set 
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of quantum Boltzmann equations describing the local particle densities and automatically 



incorporating the CP violating sources [ 125 - |127| . The closed time-path formalism is a pow- 



erful Green's function formulation for describing non-equilibrium phenomena in field theory. 
It leads to a complete non-equilibrium quantum kinetic theory approach, and a rigorous 



computation of the CP violating sources for the stop and Higgsino numbers ||1 25| |1 2Tf| . CP 
violating sources are built up when Higgsinos scatter off the advancing bubble wall, and CP 
is violated at the vertices of interactions. In the classical kinetic theory the "scattering term" 
does not include any integral over the past history of the system. This is equivalent to as- 
suming that any collision in the plasma does not depend upon the previous ones. In contrast, 
the quantum approach reveals that the CP violating source is manifestly non-Markovian. 

Let us just comment on some crucial findings. First of all, since the higgsino CP violating 
source builds up when the charginos and neutralinos scatter off the advancing bubble wall, 
there appears to be a resonant behaviour associated with the fact that the Higgs background 
is carrying a very low momentum (of order of the inverse of the bubble wall width L w ) and 
with the possibility of absorption or emission of Higgs quanta by the propagating supersym- 
metric particles ||124j| . The resonance can only take place when the higgsino and the wino 
do not differ too much in mass. By using the Uncertainty Principle, it is easy to under- 
stand that the width of this resonance is expected to be proportional to the thermalization 
rate of the -inos. A detailed computation of the thermalization rate of the supersymmetric 
particles from the imaginary part of the two-point Green function has recently been per- 
formed in ||142yi43|| by making use of improved propagators and including resummation of 



hard thermal loops, showing that the damping rate of -inos in the plasma is about 0.065 T. 
Therefore, the CP violating source for the Higgs fermion number is enhanced if M2, Mi ~ /i, 
and low momentum particles are transmitted over the distance L w . This means that the 
classical approximation is not entirely adequate to describe the quantum interference nature 
of CP violation, and only a quantum approach is suitable for the computation of how the CP 
violating sources build up. Secondly, the source is built up by integrating over all the history 
of the system. This leads to "memory" effects that are responsible for some enhancement 
of the final baryon asymmetry. These memory effects lead to "relaxation" times for the CP 
violating sources which are typically longer than the ones dictated by the thermalization 
rates of the particles in the thermal bath. 

Taking into account the detailed analysis of the bubble wall profile ||144|| , the final baryon 
asymmetry ( }4.23|) has been computed to be ||126| 



s 




4 x lO" 11 , (4.24) 



for v w ~ 1. It is intriguing that these small values of the phases are perfectly consistent 
with the constraints from the electric dipole moment of the neutron and squarks of the first 
and second generation as light as ~ 300 GeV may be tolerated. 



V. AFFLECK-DINE BARYOGENESIS 

We now turn to a third baryogenesis scenario, that has received a lot of attention. This 
mechanism was introduced by Affleck and Dine (AD) [ |145|| and involves the cosmological 
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evolution of scalar fields carrying baryonic charge. As we shall see, these scenarios are 
naturally implemented in the context of supersymmetric models. Let us quickly describe 
this before presenting a detailed analysis. 

Consider a colorless, electrically neutral combination of quark and lepton fields. In a su- 
persymmetric theory this object has a scalar superpartner, x, composed of the corresponding 
squark q and slepton I fields. 

Now, as we mentioned earlier, an important feature of supersymmetric field theories 
is the existence of "flat directions" in field space, on which the scalar potential vanishes. 
Consider the case where some component of the field x hes along a flat direction. By this we 
mean that there exist directions in the superpotential along which the relevant components 
of x can be considered as a free massless field. At the level of renormalizable terms, flat 
directions are generic, but supersymmetry breaking and nonrenormalizable operators lift the 
flat directions and sets the scale for their potential. 

During inflation it is likely that the x fi e ld is displaced from the position (x) = 0, 
establishing the initial conditions for the subsequent evolution of the field. An important 
role is played at this stage by baryon number violating operators in the potential V(x), 
which determine the initial phase of the field. When the Hubble rate becomes of the order 
of the curvature of the potential ~ rn 3 / 2 , the condensate starts oscillating around its present 
minimum. At this time, i?-violating terms in the potential are of comparable importance 
to the mass term, thereby imparting a substantial baryon number to the condensate. After 
this time, the baryon number violating operators are negligible so that, when the baryonic 
charge of x is transferred to fermions through decays, the net baryon number of the universe 
is preserved by the subsequent cosmological evolution. 

In this section we will focus on the most recent developments about AD baryogenesis, 
referring the reader to |T3| for more details about the scenario as envisaged originally in W7 



A. Supersymmetric Implementations 

The most recent implementations of the Affleck-Dine scenario have been in the context 
of the minimal supersymmetric standard model ||146| , |147|| . In models like the MSSM, with 



large numbers of fields, flat directions occur because of accidental degeneracies in field space. 
Although a flat direction is parameterized by a chiral superfield, we shall here focus on just 
the scalar component, since it is the dynamics of this that we shall be interested in. In 
general, flat directions carry a global U(l) quantum number. We shall be interested in 
those directions in the MSSM that carry B — L, since the Affleck-Dine condensate will decay 
above the weak scale and we must avoid erasure of any asymmetry through rapidly occurring 
sphaleron processes. Although there are many such directions in the MSSM ||147|| , it is 



sufficient to treat a single one as typical of the effect. As mentioned above, flat directions can 
be lifted by supersymmetry breaking and nonrenormalizable effects in the superpotential. An 
important obeservation is that inflation, by definition, breaks global supersymmetry because 
of a nonvanishing cosmological constant (the false vacuum energy density of the inflaton). In 
supergravity theories, supersymmetry breaking is transmitted by gravitational interactions 
and the supersymmetry breaking mass squared is naturally CH 2 , where C = 0(1) |148| ,P 
To illustrate this effect, consider a term in the Kahler potential of the form 
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(5.1) 



where is the field which dominates the energy density p of the universe, that is p ~ 
(/ d A 6<f)^<j)). During inflation, is identified with the inflaton field and p = V{4>) = 3H 2 m 2 . 
The term ( |5.1| ) therefore provides an effective x niass m x = —3CH 2 . This example can 
easily be generalized, and the resulting potential of the AD flat direction during inflation is 
of the form 



V{ X ) = H 2 m 2 J(x/m p ) + Hmlg( X n /< 



•pji 



(5.2) 



where / and g are some functions, and the second term is the generalized A-term coming 
from nonrenormalizable operators in the superpotential. 

Crucial in assessing the possibility of baryogenesis is the sign of the induced mass squared 
at x — 0- If the sign of the induced mass squared is negative, a large expecation value for a 
flat direction may develop, which is set by the balance with the nonrenormalizable terms in 
the superpotential. In such a case the salient features of the post-inflationary evolution are 
as follows ||147|| : 



i) Suppose that during inflation the Hubble rate is Hj ^> m 3 / 2 and the potential is of 
the form 



V( X ) = (m 3/2 - \C\Hj 
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where a is a constant of order unity, and M is some large mass scale such as the GUT 
or the Planck scale. The AD field evolves exponentially to the minimum of the potential, 
(\x\) ~ (.HjM n ~ 3 /A) 2 . The A-term in ( |5.3|) violates the U(l) carried by x, and the 
potential in the angular direction is proportional to cos(0 a + 9\ + n9) where x = |xl ei6 \ an d 
9 a and are the respective phases of a and A. This creates n discrete minima for the phase 
of x- During inflation the field quickly settles into one of these minima, establishing the 
initial phase of the field. 

%%) Subsequent to inflation, the minimum of the potential is time dependent since H 
changes with time and the AD field oscillates near the time dependent minimum 
with decreasing amplitiude. 

in) When H ~ m 3 /2 the sign of the mass squared becomes positive and the field begins 
to oscillate about x — with frequency m 3 / 2 and amplitude ~ m 3/2)- At this time 

the 5-violating A-terms are as important as the mass term, and there is no sense in which 
the baryon number is conserved. When the field begins to oscillate freely, a large fractional 
baryon number is generated during the initial spiral motion of the field. The important role 
of CP violation is also dictated by the A-terms, since the angular term cos(#a + d\ + ti9) 
becomes important at this stage, and a nonzero 9 is generated if 9 a ^ 9a- This is of course 
necessary since the baryon number is given by rif, = 2\x\ 2 9. When the condensate decays, 
the baryon asymmetry is transferred to fermions. 

iv) In the case in which the inflaton decays when H < m 3 / 2 (consistently with the 
requirement that not too many gravitinos are produced at reheating), the resulting baryon 
to entropy ratio is 
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and T r is the reheating temperature. For example n = 4 gives ||147 



(5.4) 



where n x = p x /m 3 / 2 , p^ ~ m l/2 m l is the energy density of the inflaton field when H ~ rn, 3 / 2 , 



B ~ 10 
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10 6 GeV, 



10~ 3 M 



(5.5) 



The MSSM contains many combinations of fields for constructing flat directions carrying 
a nonvanshing B — L [ 147|| . Particularly appealing directions are the ones which carry B — L 
and can be lifted at n = 4 level. They are the x 2 = LH2 directions. The nonrenormalizable 
operator is then 



A 
M 



W = -- (LH 2 ) 2 



(5.6) 



which may be present directly at the Planck scale, or could be generated, as in 50(10) 
GUTs, by integrating out right-handed neutrino superfields which are heavy standard model 
singlets. 



B. Initial conditions for viable AD baryogenesis 

^From the above picture it is clear that a successful AD baryogenesis mechanism is 
achieved if the effective potential during inflation contains a negative effective mass term 
and nonrenormalizable terms that lift the flat directions of the potential. 

With a minimal Kahler potential only, the effective mass squared m 2 ~ H 2 during infla- 
tion is positive, x = is stable, and the large expectation values required for baryogenesis 
do not result. Quantum de Sitter fluctuations do excite the field with (x 2 ) ~ H], but with 
a correlation length of the order of Hj . Any resulting baryon asymmetry then averages 
to zero over the present universe. In addition, the relative magnitude of the B violating 
A-terms in the potential is small for Hi <^ M. 

A negative sign for the effective mass squared is possible if one either considers general 
Kahler potentials, or in the case in which inflation is driven by a Fayet-Iliopoulos D term, 
which preserves the flat directions of global supersymmetry, and in particular keeps the 
inflaton potential flat, making D-term inflation particularly attractive ||79|| . If the AD field 



carries a charge of the appropriate sign under the U(l) whose Fayet-Iliopoulos term is 
responsible for inflation, a nonvanishing vacuum expectation value may develop |[L50|| . Notice 
that, if the AD field is neutral under the U(l), it may nonetheless acquire a negative mass 
squared after the end of D-term inflation through nonminimal Kahler potentials relating 
the AD field to the fields involved in the inflationary scenario | |151| . 



Another (marginal) possibility is that the effective mass squared during inflation is pos- 
itive, but much smaller than H], < C < 1. This happens either in supergravity models 
that possess a Heisenberg symmetry in which supersymmetry breaking by the inflationary 



vacuum energy does not lift flat directions at tree level [|1 52| , or in models of hybrid inflation 



based on orbifold constructions, in which a modulus field T is responsible for the large value 
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of the potential during inflation, and a second field with appropriate modular weight is re- 
sponsible for the roll-over [|150|| . The correlation length for de Sitter fluctuations in this case 
is i c = Hi exp(3H]/2m 2 ). This is large compared to the horizon size. In fact, the present 
length corresponding to i c should be larger than the horizon size today. Using Hi = 10 13 
GeV, this gives (m/Hi) 2 <> 40, which is easily satisfied. The flat direction is truly flat and 
the AD baryogenesis may be implemented as originally envisaged f|7j . After inflation, these 
truly flat directions generate a large baryon asymmetry, typically B = 0(1). Mechanisms 
for suppressing this asymmetry to the observed level have been considered in ||152| . These 



include dilution from inflaton or moduli decay, and higher dimensional operators in both 
GUT models and the MSSM. The observed BAU can easily be generated when one or more 
of these effects is present. 



C. Preheating and AD baryogenesis 

As we have described in subsection ( |1UB| ), in the first stage of reheating after inflation, 
called preheating, resonance effects may lead to explosive particle production. Fluctuations 
of scalar fields produced at the stage of preheating are so large that they may break super- 
symmetry much more strongly than inflation itself. These fluctuations may lead to symmetry 
restoration along flat directions of the effective potential even in the theories where the usual 
high temperature corrections are exponentially suppressed. These nonthermal phase tran- 
sitions after preheating can play a crucial role in the generation of the primordial baryon 
asymmetry by the AD mechanism. In particular, the baryon asymmetry may be generated 
at the very early stage of the evolution of the universe, at the preheating era, and not when 
the Hubble parameter becomes of order the gravitino mass 7713/2 ||153| . 



Let us consider the chaotic inflation scenario where the inflaton field couples to a 
complex X-field with a potential of the form 

M 2 

V = ^<j ) 2 + g 2 <j ) 2 \X\ 2 . (5.7) 

After inflation, the initial energy of the inflaton field is transferred to the X-quanta in the 
regime of parametric resonance |87| . It is especially important that the maximum amplitude 



of field fluctuations produced at that stage may be sizeable, (X 2 ) ~ M^/g 2 or even 
larger. These large fluctuations occur only in the bosonic sector of the theory, thus breaking 
supersymmetry by a large amount at the quantum level. Suppose the field x labelling the 
AD flat direction couples to the generic super mult iplet A by a renormalizable interaction 
of the form £ int = h 2 \x\ 2 \X\ 2 . Preheating gives an additional contribution to the effective 
squared mass along the x-direction beyond the contribution of the — \C\H 2 term discussed in 
the previous subsection, which is responsible for a nonvanishing vacuum expectation value 
(Xo) of the AD field during inflation. The curvature of the effective potential becomes large 
and positive, and the symmetry rapidly restores for h 2 (X 2 ) > \C\ H 2 . Thus, for C = 0(1), 
symmetry becomes restored for h 2 > g 2 . In other words, supersymmetry breaking due to 
preheating is much greater than the supersymmetry breaking proportional to H unless the 
coupling constant h 2 is anomalously small. This leads to symmetry restoration in the AD 
model soon after preheating. The AD field rapidly rolls towards the origin x — and makes 
fast oscillations about it with an initial amplitude given by x ~ Ixol- The frequency of 
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the oscillations is of order of h(X 2 ) 1 ^ 2 and is much higher than the Hubble parameter at 
preheating. The field is underdamped and relaxes to the origin after a few oscillations. 
Baryon number can be efficiently produced during this fast relaxation of the field to the 
origin. This picture is similar to that of ref. | i7fl , but there are considerable differences. 



According to |47[], the baryon asymmetry is produced at very late times when the Hubble 
parameter H becomes of order of the gravitino mass m.3/2. At preheating, baryon number 
production occurs right after inflation. This is possible because the same nonthermal effects 
generating a large curvature at the origin give rise also to sizable CP- violating terms along 
the AD flat direction 



D. Q-balls and AD baryogenesis 

Q-balls are stable, charged solitons appearing in a scalar field theory with a spontaneously 



broken global U(l) symmetry | 154 |. The Q-ball solution arises if the scalar potential V((j)) is 
such that V (<p) / \<f)\ 2 has a minimum at non-zero 0. In the MSSM such solitons are part of the 
spectrum of the theory [ |155|| . This is due to the form of the MSSM scalar potential and the 
fact that squarks and sleptons carry a global B — L charge. An interesting possibility for Q- 



ball formation is provided by the fact that there are many flat directions in the MSSM | 1 17| . 



During inflation the MSSM scalar fields are free to fluctuate along these flat directions and 
to form condensates. This is closely related to Affleck-Dine (AD) baryogenesis [H7J with the 
difference that the condensate is naturally unstable with respect to the formation of Q-balls 
|, |155| with a very large charge yi57[| . The properties of the MSSM Q-balls depend 



upon the scalar potential associated with the condensate, which in turn depends upon the 
supersymmetry breaking mechanism and on the order at which the non-renormalizable terms 
lift the degeneracy of the potential. If supersymmetry breaking occurs via the supergravity 



hidden sector ||45|| , B-balls are allowed to form if the charge density inside the initial lump 
is small enough. Since this charge density scales as the inverse of the reheating temperature 
T r , this can be translated into the condition [|158| , |159| 



On the other hand, after formation, B-balls could be dissociated by the bombardment of 
thermal particles, or dissolve by charge escaping from the outer layers. Both phenomena 
may be prevented if T r <^ (10 3 — 10 5 ) GeV. Therefore, surviving B-balls must have very large 
charges, B <; 10 14 ||158|| . Since the decay rate of the B-ball also depends on its charge, for 
B <; 10 14 , B-balls will indeed decay well below the electroweak phase transition temperature. 
The only requirement is relatively low reheating temperature after inflation, typically of the 
order of 1 GeV, which is fixed by the observed baryon asymmetry when the CP violating 



phase responsible for the baryon asymmetry is of the order of unity |159| . In particular, this 



is expected to be true for D-term inflation models ||151| , |160|| . Such B-balls can therefore decay 



at temperatures less than that of the electroweak phase transition, and consequently lead to a 
novel version of Affleck-Dine baryogenesis, in which the baryon asymmtry asymmetry comes 
from B-ball decay rather than condensate decay. Since the baryon asymmetry production 
takes place when sphaleron processes are no longer active, this mechanism can work even in 
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the presence of additional L violating interactions or B — L conservation, which would rule 
out conventional AD baryogenesis along B violating directions. 

VI. CONCLUSIONS 

The origin of the baryon to entropy ratio of the universe is one of the fundamental initial 
condition challenges of the standard model of cosmology. Over the past several decades, 
many mechanisms utilizing quantum field theories in the background of the expanding uni- 
verse have been proposed as explanations for the observed value of the asymmetry. In this 
review, with space constraints, we have focused on those that we feel hold the most promise. 
We have given what we hope is an up to date review of three popular theories of baryogenesis; 
GUT scenarios, electroweak baryogenesis, and the Affleck-Dine mechanism. Each of these 
mechanisms has both attractive and problematic aspects, as we'll attempt to summarize 
here. 

Grand Unified baryogenesis was the first implementation of Sakharov's baryon number 
generation ideas. There are at least three reasons why this is an attractive mechanism. 
First, baryon number violation is an unavoidable consequence of Grand Unification. If 
one wishes to have quarks and leptons transforming in the same representation of a single 
group, then there necessarily exist processes in which quarks convert to leptons and vice- 
versa. Hence, baryon number violation is manifest in these theories. Second, it is a simple 
and natural task to incorporate sufficient amounts of CP-violation into GUTs. Among other 
possibilities, there exist many possible mixing phases that make this possible. Finally, at 
such early epochs of cosmic evolution, the relevant processes are naturally out of equilibrium 
since the relevant timescales are slow compared to the expansion rate of the universe. The 
simplicity with which the Sakharov criteria are achieved in GUTs makes this mechanism 
very attractive. 

While GUT baryogenesis is attractive, it is not likely that the physics involved will be 
directly testable in the foreseeable future. While we may gain indirect evidence of grand 
unification with a particular gauge group, direct confirmation in colliders seems unrealistic. 
A second problem with GUT scenarios is the issue of erasure of the asymmetry. As we 
have seen, the electroweak theory ensures that there are copious baryon number violating 
processes between the GUT and electroweak scales. These sphaleron events violate B + L 
but conserve B — L. Thus, unless a GUT mechanism generates an excess B — L, any baryonic 
asymetry produced will be equilibrated to zero by anomalous electroweak interactions. While 
this does not invalidate GUT scenarios, it is a constraint. For example, SU(5) will not be 
suitable for baryogenesis for this reason, while SO(10) may be. This makes the idea of 
baryogenesis through leptogenesis particularly attractive. 

In recent years, perhaps the most widely studied scenario for generating the baryon num- 
ber of the universe has been electroweak baryogenesis. The popularity of this idea is tightly 
bound to its testability. The physics involved is all testable in principle at realistic colliders. 
Furthermore, the small extensions of the model involved to make baryogenesis successful can 
be found in super symmetry, which is an independently attractive idea, although electroweak 
baryogenesis does not depend on supersymmetry. 

However, the testability of electroweak scenarios also leads to tight constraints. At 
present there exists only a small window of parameter space in extensions of the electroweak 
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theory in which baryogenesis is viable. This is because electroweak baryogenesis to be 
effective requires a strong enough first order phase transition. This translates into a severe 
upper bound on the lightest Higgs boson mass, m h <> 105 GeV, in the case in which the 
mechanism is implemented in the MSSM. Furthermore, the stop mass may be close to 
the present experimental bound and must be smaller than, or of order of, the top quark 
mass. The most direct experimental way of testing this scenario is through the search for 
the lightest Higgs at LEP2. In this sense, we are close to knowing whether electroweak 
processes were responsible for the BAU. 

If the Higgs is found, the second test will come from the search for the lightest stop at 
the Tevatron collider (the stop mass is typically too large for the stop to be seen at LEP). 
If both particles are found, the last crucial test will come from B physics, more specifically, 
in relation to the CP-violating effects. Depending on their gaugino-Higgsino composition, 
light neutralinos or charginos, at the reach of LEP2, although not essential, would further 
restrict the allowed parameter space consistent with electroweak baryogenesis. 

Moreover, the selected parameter space leads to values of the branching ratio BR(5 — > 57) 
different from the Standard Model case. Although the exact value of this branching ratio 
depends strongly on the value of the fi and A t parameters, the typical difference with respect 
to the Standard Model prediction is of the order of the present experimental sensitivity and 
hence in principle testable in the near future. Indeed, for the typical spectrum considered 
here, due to the light charged Higgs, the branching ratio BR(6 — > S7) is somewhat higher 
than in the SM case, unless negative values of A t ji are present. 

Affleck-Dine baryogenesis is a particularly attractive scenario, and much progress has 
been made in understanding how this mechanism works. As was the case for electroweak 
baryogenesis, this scenario has found its most promising implementations in supersymmetric 
models, in which the necessary flat directions are abundant. Particularly attractive is the 
fact that these moduli, carrying the correct quantum numbers, are present even in the 
MSSM. 

The challenges faced by Affleck-Dine models are combinations of those faced by the GUT 
and electroweak ideas. In particular, it is necessary that B — Lbe violated along the relevant 
directions (except perhaps in the Q-ball implementations) and that there exist new physics 
at scales above the electroweak. If supersymmetry is not found, then it is hard to imagine 
how the appropriate flat directions can exist in the low energy models. 

Understanding how the baryon asymmetry of the universe originated is one of the fun- 
damental goals of modern cosmology. Particle physics has provided us with a number of 
possibilities, all of which involve fascinating physics, but which have varying degrees of 
testability. Our primary goal in this review has been to provide an as up to date as possible 
account of current efforts in three primary avenues of baryogenesis research. We have also 
tried to be as pedagogical as space allows. We hope that the reader has found this helpful, 
and that friends and colleagues will forgive us for any inadvertent omissions. 
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